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Abstract 



Stationary wave solutions of the perturbed Korteweg-de Vries equation are considered in 
the presence of external hamiltonian perturbations. Conditions of their chaotic behaviour 
are studied with the help of Melnikov theory. For the homoclinic chaos Poincare sections are 
C^l ■ constructed to demonstrate the complicated behaviour, and Lyapunov exponents are also 

numerically calculated. 
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G\ ■ 1 Introduction 

o , 

Recently there has been a significant attention to the chaotic behaviour of the solutions to partial 
differential equations. For example, chaos was found in a perturbed sine-Gordon equation 
nonlinear Schrodinger equation Q , and later for Korteweg-de Vries-Burgers equation Q and its 
generalization on high-order nonlinearities and the case of Kadomtsev-Petviashvili equation y] . 
It was shown how chaos can appear in such systems which were completely integrable without 
perturbation via the appearance of subharmonics and homoclinic tangles. For the KdV equation 
incorporating dissipation and instability || it was shown that for the strongly dissipative case 
the overall evolution of solutions is chaotic with irregular soliton interactions. 

In our previous paper Q we've considered the influence of time-periodic hamiltonian external 
perturbations on a KdV system and showed, how the stochastic layer can appear on a phase 
plane near the unperturbed separatrix. The width of this layer was calculated with the help 
of Chirikov criterion for the overlap of resonances. In this work we continue this study on 
an example of one-harmonic perturbation. The fact that this perturbation is taken explicitly 
provides the possibility to obtain an exact expression for the Melnikov function in order to 
determine the transition to a chaotic behaviour. To study the system in a near-separatrix 
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Figure 1: (Left) Phase plane of the system (|2.1|). (Right) Plot of the potential energy ( |2.2| ). 
Here v = 1 and (3 = 1. 

region, we transform the governing ODE into a map for which the width of stochastic layer can 
be easily obtained. Finally we plot Poincare sections of chaotic behaviour and calculate the 
corresponding Lyapunov exponents. 

We start with the perturbed Korteweg-de Vries equation taken in the form: 

ut + cu x + uu x + (3u xxx = f(u, x - Vt) x . (1.1) 

Here f(u, x — Vt) is assumed to be periodic on its last argument and will be taken simply as 
f = fo coslo(x — Vt). Transforming coordinates in a moving frame (x 1 — > x — Vt, t' — ► t) and 
considering steady waves Ut = 0, one obtains (the primes are omitted): 

u 2 

fiu xx = vu - — + f cos ojx + C, (1.2) 

where V = c + v (supercritical case) or V = c — v ( subcritical case), v > 0, C is the integration 
constant which can be chosen equal to zero from the condition that the solution is bounded at 
infinity. Without loss of generality we consider here only the supercritical case. Therefore we 
have: 

u 2 

f3uxx = vu - — + fo coswx (1.3) 

or with z = u x 

f 

u x = z 

M \ (L4) 

Zx = p \ vu ~ \ + fo cos UJX ) ■ 
The outline of the rest of the paper is as follows. In the next Section the unperturbed case of 

the Korteweg-de Vries equation is considered and the explicit expressions for the solutions are 

obtained. Later in Section 3 conditions of chaos are derived on the basis of studying subharmonic 



and homoclinic Melnikov functions. Then, in Section 4, we transform our ODE ( |1.3| ) into a 
mapping, for which phase plots are presented together with the calculation of the width of 
the stochastic layer. Section 5 contains Poincare plots numerically found in the vicinity of a 
saddle along with the Lyapunov exponents which prove to be positive thereby supporting our 
conclusion about chaotic behaviour in the system. Finally, Section 6 contains a summary and 
conclusions. 



2 Unperturbed case 



For the unperturbed system (fo = 0) we simply have 



/3u xx = vu 



u~ 



(2.1) 



This is an equation of motion for a nonlinear oscillator in a potential field 
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with the Hamiltonian 



P( ") = T-- 



TT . . (3z 2 u 3 vu 2 
H {u, z, x) = — + — - — . 



(2.2) 



(2.3) 



If we now look for the solutions of (2.1) in the form of elliptic Jacobi functions of modulus to: 



u m (x) = a + bcn 2 (kx \ to), < to < 1, 



(2.4) 



where a, b and k are depending on to constants, then the following expressions can be easily 
obtained: 



u m (x) = v+ < x ~ 2m \ + 

z m (x) = Svrnk 

^ 1 %/m 2 — m+1 

T m _ 2K(m) 
1 ~ k ' 



3vm cn 2 {kx | to), 



+1 \J m?—m+l 

cn.(kx | m)sn(kx \ to) dn(kx \ to), 



(2.5) 



k = (to 2 -to + 

Here T m is the period of oscillations and K{m) is the complete elliptic integral of the first kind. 
In the limit m — ► 1 fl2.5|) yields the separatrix 



Uq(x) = 3v sech 2 (|y / ^'a;), 

zq(x) = — 3uw ^ sinh(^. /^x) sech 3 (^. /^x). 



(2.6) 



Phase plot of the system fl2,l|) as well as the potential energy (|2,2|) is presented in Figure 1. 
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Figure 2: The width of stochastic layer as the function of perturbation frequency Q4.5|). Here 



P = 1, v = 1, / = 0.01. 



3 Melnikov chaos 



The source of chaotic motion in our system is the homoclinic orbit on a phase plane. In the 
unperturbed case this orbit is formed by coinciding stable and unstable manifolds of the sad- 
dle. When there is a perturbation, the homoclinic orbit can be broken to yield the transverse 
intersection of stable and unstable manifolds, which gives rise to chaotic behaviour near the 
separatrix. As far as one of the precursors of chaos is the appearance of subharmonics, we'll 
start with the subharmonic Melnikov function defined for the periodic orbits as 



with the period T m = rT/s for r, s relatively prime integers, g is the small vector of perturbation, 
and the wedge product is defined by a A b = ai&2 — a^i- 

In considering our problem we substitute the term z m (x) for /, and the term of perturbation 
/o cos ujx for g into the subharmonic Melnikov function, therefore obtaining: 





where f(q m ) is the unperturbed vector field corresponding to the solution q m (x) 



(u m ,z m ) T 





Calculations with the help of Fourier representation of elliptic functions 



finally give: 




(3.3) 



where 




(3.4) 
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Figure 3: Separatrix mapping (gj) with v = 1, (3 = 1, u = 3. (Left) / = 1.83 • KT 3 . (Right) 
/o = 9.1 • 1(T 3 . 



As it is clearly seen from M r / S {6) has simple zeros, what proves the presence of subharmonic 
bifurcations ||. 

To move further the homoclinic Melnikov function M{6) can be introduced as a limit of 
M r / S (9) when s = 1, r —* oo, and m — ► 1: 

M(0) = ^ /°° *b(aO cos [w(x + 0)] da; = 12mj2 -^_ sin ^. (3.5) 
^ sinh7ru;\/5 



Surely this function also has simple zeros on 9 and therefore proves the existence of transverse 
homoclinic orbits. The presence of these orbits means, via Smale-Birkhoff theorem, the ap- 
pearance of a Smale horseshoe in the vicinity of an unperturbed saddle and thereby homoclinic 
chaos in this region of a phase space ||. We note here that unlike the situation with dissipation 
( like in || or Q), when Melnikov theory resulted in some relation between the perturbation 
amplitude, frequency and the dissipation coefficient, in our case there are no such restrictions. 
It should be also noted that although the presence of a stochastic layer near the separatrix can 
be stated for rather generic hamiltonian perturbation making system near-integrable [||, the 
transverse intersection of stable and unstable manifolds is still a question which should solved 
by Melnikov theory in every particular case. 



4 Near-separatrix motion 

To study the chaotic motion near the separatrix we shall use the whisker mapping (separatrix 
mapping) |]l0| . To derive it we note first that the change in the unperturbed energy is given by: 

dH dH dV 



dt dp dx 

5 



Pfoz(x) cos ujx, (4.1) 
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Figure 4: Separatrix mapping 
f = 9.1 • 1(T 2 . 



with v = 1, /3 = 1, U3 = 3. (Left) / = 1.83 • 1(T 2 . (Right) 



where p = (3z is the momentum. Since the perturbation is a periodic function of time, we may 
introduce the phase of perturbation: 



ujx + const, 



(4.2) 



which is an equivalent of time. In order to obtain the desired separatrix mapping one has to 
consider the discretized time scale x n , and the mapping will involve <j) n and E n as variables. 
The change of phase after one period of oscillations is equal to LvT m (E n+ i), while energy change 
AE n+ \ = E n + AE n . Now we can estimate the change in the unperturbed energy AE n per 
period of motion in proximity separatrix at time x n : 

r x„+T/2 



AE n = Pf I z(x) cos [oj(x + x n )} dx « /3f 

Jx n -T/2 



zq(x) cos [u(x + x n )\ dx. 



(4.3) 



Here we have approximated AE n by evaluating the integral (4.3) on the unperturbed separatrix 
in accordance with the standard procedure 0]. As it is clearly seen from the comparison of 
(|3.5|) and (|4.3|), AE n = (3 2 M(x n ). Expansions of T m and E(m) in the vicinity of m = 1 give 



finally the following correlation between them: T(E) as 2^J (v/(3) In 24/^/2 | E \ . Using this 
approximate period we can write the separatrix mapping in the following form: 

12 ^ 2 r sin^, 

(4.4) 



E n +i 

4>n+l 



E n + 



sinh i 



> sin i 

24 



6 « + 2 v / ^ ln 7fe (mod27r) - 



If to measure the stretching of a small phase interval with the parameter K = \8(f) n +i / 5(f) n — 1| 
§, one can detect the border of chaotic motion as an appearance of a local instability in phase: 
K > 1. For the mapping (^4.4|) this condition gives the following approximation for the width of 
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Figure 5: (Left) Poincare section of the system (p..3|). (Right) The same near the saddle (0,0). 
Here fo = 0.1, u = 15, v = 1, /3 = 1, and initial conditions are taken in the point (0.0001,0). 
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Figure 6: The same as in Figure 5 with to = 20 



the stochastic layer E s t- 



rwi2irf uj 3 p 2 
ft sinh ttuj ,/-| 



E\< E 



St 



(4.5) 



Typical form of this dependence is presented in the Fig. 2. In Fig. 3, 4 we have plotted the 
phase plane of the mapping (E, <j>) for different values of the perturbation amplitude /o- Fig. 
3(Right) shows how the stochastic layer arises around the unperturbed separatrix. Then the 
width of this layer increases together with fo, leading finally to a chaotic sea, as demonstrated 
in Fig. 4(Right). 



5 Numerical simulations 



In this section we present the Poincare sections of (F3) in the vicinity of the saddle (0, 0), which 
show homoclinic crossings in this region. The methods used were Runge-Kutta-Fehlberg and 
Adams-Bashford-Moulton ("predictor-corrector") method. 
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Figure 7: The same as in Figure 5 with u> = 25 



Local chaotic behaviour is detected by the numerical calculation of the dominant Lyapunov 



exponent. It was computed using the method due to [11]. The system of equations, which 
defines the value of this exponent, was integrated with the Runge-Kutta fourth-order algorithm, 
while the equation fll.3| ) for a fiducial trajectory was run with fourth-order symplectic scheme 
[12]. In Figures 5-8 the Poincare sections are presented for the system ( |1.3| ) with the forcing 
amplitude /o = 0.1 and different values of the frequency: uo = 15, 20, 25 and 30. The value of 
the largest Lyapunov exponent in these cases is A = 0.107, A = 0.149, A = 0.169 and A = 0.188 
respectively, and so we can state that we really have chaotic motion in a close vicinity of a 
separatrix, as it was predicted. Numerical calculation of the Lyapunov exponent in the region 
of regular dynamics with initial conditions at (1.5, 0) yields 0.0001, and therefore our results can 
be taken as significant. 



6 Conclusions 

In this work we have considered the stationary wave reduction of the Korteweg-de Vries equation 
under small hamiltonian perturbations. On the basis of Melnikov theory it was shown how 
chaos can occur in this situation via the appearance of subharmonics and a further homoclinic 
tangle. Theoretical predictions about chaotic behaviour are supported numerically by plotting 
Poincare sections and calculation of the corresponding Lyapunov exponents. In a near-separatrix 
region the governing ODE was transformed into a mapping what allowed one to find the width 
of the stochastic layer. 

Our future work will focus on the development of kinetic description of chaotic behaviour for 
the obtained mapping || , and with the derivation of new methods for the symplectic integration 
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Figure 8: The same as in Figure 5 with u> = 30 



of the KdV equation under hamiltonian perturbations following the framework of [13], where 
the same was done for sin-Gordon system. 
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